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ABSTRACT 



The conceptual development of a computer flight simulation for 
design, testing and analysis of departure prevention systems, simu- 
lation capability and programming are discussed, along with required 
research material and data. A description is given of the aerodynamic 
buildup program written for incorporation in the simulation, including 
the aerodynamic equations of the model base aircraft, sample program 
statements and output. 
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I. INTRODUCTION 



Throughout the history of aviation, departure from controlled flight 
has been a persistent problem. Departure has occurred during various 
periods of aviation history for different reasons. In the early years, 
it was an inadequate knowledge of aerodynamic effects leading to poor 
or inadequate designs. In more recent years, modern design techniques 
and an improved understanding of aerodynamics, and stability and control 
have led to the design of high performance aircraft which constantly fly 
at the limits of their operating envelopes and that in less than a 
seconds time can be outside of that envelope departing controlled flight. 
In past times, recovery from departure was often a relatively easy pro- 
cedure. It still is with simple, basic, fundamental, stable aircraft 
designs. Recent state-of-the-art tactical aircraft, however, realize 
their capabilities by displaying neutral or unstable static stability 
compensated for by digital fly-by-wire control systems. These aircraft 
with their instabilities and non-conventional aerodynamic design features 
are not so easily recovered. 

As aircraft control systems have been developed over the years, many 
and varied departure control , departure prevention and departure recovery 
systems have been developed and flown. The majority of these systems 
have been limiting type systems, which in some way limit the operation 
of the aircraft; an angle-of-attack limiter being a common example. 

During the performance of an aircraft mission, an actual departure, 
whether controlled or uncontrolled, recoverable or unrecoverable, will 
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result in at least the loss of mission effectiveness and probably the 
loss of man or aircraft or both. By the same means, restructing air- 
craft operation to levels below the maximum designed capability in order 
to avoid potential departure situations may result in the same losses of 
mission, man and/or aircraft. For these reasons it is desirable to de- 
velop a departure prevention system for tactical aircraft that is as 
"non-limiting" as possible. 

This thesis is the first report on the development of a computer 
flight simulation for the design, testing and analysis of modern optimal, 
adaptive departure systems. It contains the results of project defini- 
tion and planning, and the details of the aerodynamic buildup developed for 
incorporation in the flight simulation program package. 
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II. SIMULATION CONCEPT DESCRIPTION 



A. SIMULATION CAPABILITY 

The development of a flight simulation is very dependent on the 
purpose for which it will be utilized. A full flight simulation is 
required for full motion base simulator, whereas a much reduced version 
may be used for investigation of carrier landing characteristics. The 
following are some of the key points considered and decisions made in 
determining the type and extent of the simulation needed for this 
project. 

1. Although data indicates that departure is still a problem in 
older tactical aircraft, the application of modern active control tech- 
niques to departure systems is most applicable to fly-by-wire or control - 
by-wire systems. 

2. Availability of data led to utilization of the McDonnell Douglas 
F/A-18A as the simulation data base. 

3. The desire to avoid the additional knowns and unknowns of super- 
sonic flight performance reduced the simulation speed envelope to the 
subsonic regime. 

4. For the most applicable case, the simulation will involve up- 
and-away flight conditions only. 

5. The outer loop closures of the aircraft automatic flight control 
system will not be simulated but in its place an outer-loop maneuvering 
autopilot will be modeled. The aircraft control augmentation system 
will be simulated. 
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6. Given the above conditions and the potential to depart flight 
throughout the entire flight envelope the full aircraft system in terms 
of operating limits, control laws and systems will be modeled as closely 
as possible to the model base aircraft. 

The resulting flight simulation will be comparable with other digi- 
tal fly-by-wire aircraft. Controlled maneuvers will be precisely per- 
formed and repeatable via the maneuvering autopilot and the performance 
and flying qualities should match closely with that of the F/A-18. 

B. SIMULATION PROGRAMMING 

The programming of a flight simulation generally consists of three 
major components, a flight control laws model, an aerodynamic buildup, 
and flight dynamics calculations. Each of these components is quite 
complex in itself with the entire simulation requiring several program- 
mers. This results in a modular type programming with each of the three 
components comprising a module. This is an optimum situation in that 
each module, control laws, aerodynamic buildup, and flight dynamic per- 
forms different calculations for which programming can be specifically 
tailored. Once programmed, each module can be tested by test stubs to 
verify results prior to inclusion in the full flight simulation program. 
The use of modular programming reduces the complexity of the simulation 
and allows identification of real and potential problems in the simula- 
tion by testing each module through the full range of flight conditions. 
The tailoring of the programming for the various modules led to the 
utilization of both CSMP and FORTRAN computer languages, in the simula- 
tion. The appropriate language is utilized in the simulation where the 
following characteristics are advantageous; 
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CSMP: 



- The capability to handle nonlinear and time-invariant problems. 

- The provisions to allow the modeling/simulation of a physical 
system utilizing block diagrams. 

FORTRAN : 

- The capability to handle a large quantity of data. 

- The capability for formatted output. 

- The capability for logic, branching and subroutines. 

CSMP is generally used for the dynamic flight control laws model, the 
flight dynamics calculations, the program controls and unformatted out- 
put. FORTRAN language is used for the aerodynamic buildup, gain func- 
tions, other minor functions where necessary and the simulation formatted 
output. 

The use of each language where appropriate results in a faster, more 
accurate, more efficient flight simulation. 

C. SIMULATION FORMAT AND OPERATION 

The flight simulation program format consists of the three major 
modules; flight control laws model, aerodynamic buildup and flight 
dynamics calculations along with a program explanations section, a pro- 
gram control and output section and minor subroutines and functions. 

The following is a brief description of the operation of the flight 
simulations program (see Figure 1). Command inputs are made to the 
dynamic flight control laws model (CSMP), the output of which are con- 
trol surface defections, of leading edge and trailing edge flaps, ailer- 
ons, horizontal stabilizer, rudder and speedbrake. These surface deflec- 
tions are input to the aerodynamic buildup (FORTRAN) the output of which 
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are the aircraft total coefficients for lift, drag, pitching moment, 
rolling moment, yawing moment and side force. These coefficients are 
inputs to the flight dynamics module where aircraft rotational and trans- 
lational motions of pitch rate, roll rate, yaw rate and U, V, W velocities 
are computed from the equations of motion. The aircraft motions are fed 
back to the command input side of the flight control laws module for com- 
parison to commanded inputs and subsequent command modification. The 
program run time, integration time and other control functions are input 
from the CSMP program. Output is generated from both CSMP statements 
and FORTRAN subroutines for formatting. 
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III. DATA AND RESOURCES 



In developing a flight simulation, two types of information are 
needed: (1) required information - flight systems description, etc., 

and (2) reference information - programming options, etc. Reviewing 
wide range of tasks required for a simulation of this magnitude the 
need to have a source library is obvious. The project has four dis- 
tinct tasks to be performed, (1) project definition and planning, (2) 
mathematical modeling, (3) programming and (4) testing and analysis. Re- 
search material and required data was collected in each of these areas 
for use in completing the tasks. The collected material can be divided 
into six areas, (1) General Departure Information, (2) Aerodynamic Data, 
(3) Flight Control Laws, (4) Maneuvering Autopilots, (5) Programming 
Techniques, and (6) Flying Qualities. The following is a list of the 
major resources obtained for the flight simulation project, and a brief 
description of each. 

A. GENERAL DEPARTURE INFORMATION 

Reference 1 contains all mishap reports from mishaps classified by 
type as uncontrolled flight. It is subdivided into jet, prop and heli- 
copter mishaps and provides information on mishap causes, phase of mis- 
sion and a narrative of the mishap. 

B. FLIGHT CONTROL LAWS 

Reference 2 is a description of the inner and outer loop control 
laws. It is presented in three sections as follows: 
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1. Flight Control System Characteristics: Inner Loop Theory of 

Operation. This section contained information on the longitudinal, 
lateral and directional control laws, quad sensor signals, actuator 
systems, angle-of attack system and air data system. 

2. Automatic Flight Control System: Theory of Operation. 

3. Autothrottle: Theory of Operation. 

Reference 3 contains system descriptions and diagrams of the follow- 
ing systems pertinent to a flight simulation: longitudinal and lateral- 

directional control systems, flap commands, mechanical primary controls, 
flight control electronic set, actuation devices, and throttle control. 

Reference 4 is the F/A-18 version 8.2.1 flight control system descrip- 
tion and theory of operation. It contains a description of the flight 
control hardware and interfaces and the system theory of operation in- 
cluding software architecture and mathematical characteristics of inner 
and outer loop control laws. 

C. STABILITY AND CONTROL 

Reference 5 contains the stability and control characteristics of 
the production F/A-18 high speed maneuvering and high lift configurations, 
derived from wind tunnel test and revised where appropriate to reflect the 
results of developmental flight tests. The report presents data in a 
graphical form for status longitudinal and lateral-directional stability 
and control, and the longitudinal, and lateral -directional dynamic 
derivatives. 

D. FLYING/HANDLING QUALITIES 

Reference 6 presents the flying and handling qualities of the F/A-18 
figher escort configuration. Longitudinal and lateral-directional modes 
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and responses, unaugmented characteristics, and spin departure character- 
istics are included. Information is in both graphical and tabular form. 

E. MANEUVERING AUTOPILOT 

Reference 7 is a discussion of developing a maneuvering autopilot. 

It includes maneuvering requirements, linear analysis and design, control 
law development, command generation and flight experience. 

F. PROGRAMMING TECHNIQUES 

Information on programming techniques for manipulating large qualities 
of data with emphasis on flight simulations and aerodynamic buildups was 
obtained from both Northrop Aircraft Corp. and the Naval Air Development 
Center. 

This is by no means a complete list of the information obtained. It 
is, however, the primary material used during the project. It is dis- 
cussed to indicate the type of materials required to develop a flight 
simulation. The general departure material was used to determine what 
flight conditions should be investigated. The flight control laws mater- 
ial is being utilized to develop the dynamic flight control law model. 

The stability and control data is used in the aerodynamic buildup. The 
maneuvering autopilot data is used for modeling the outer loop maneuver- 
ing autopilot. The programming techniques material is used for program- 
ming methodology and the flying qualities data is used for verification 
of simulation model response. 
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IV. AERODYNAMIC BUILDUP 



A. CONSIDERATIONS 

As discussed earlier, the major parts to a flight simulation program 
are a flight control laws model, an aerodynamic buildup and flight 
dynamics calculations. The following is a description of the aerodynamic 
buildup developed for incorporation into the flight simulation program. 

In developing the buildup, the following goals were set. 

1. Simplicity and intelligibility. 

2. Ability to operate as a separate program or be incorporated as 
a subprogram in a larger simulation. 

3. Provide proper results throughout the entire range of flight con- 
ditions for the simulation. 

4. Flexibility, versatility and alterabil ity. 

The aerodynamic buildup constitutes a large portion of the entire simula- 
tion. It also involves the manipulation of very large quantities of data. 
Its programming must consider, integration with other program modules, 
data handling times and storage space. These considerations impact on 
decisions about programming language, programming methodology, and data 
storage and retrieval techniques. 

B. AERODYNAMIC EQUATIONS 

The operation of the flight simulation program, discussed in Chapter 
Two, indicated the inputs to the aerodynamic buildup are the aircraft 
control surface deflections and the outputs are the aircraft aerodynamic 
coefficients. The first task was the determination of what control 
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surface deflections and flight conditions affected each coefficients and 
to what extent. For example, lift coefficient is changed by deflecting 
the horizontal stabilizer. How much it is changed is determined by the 
amount of deflection, the airspeed, and the angle-of-attack. This infor- 
mation was determined from the model base aerodynamic equations [Ref. 4], 
Below is a list of the control surfaces and flight conditions affecting 
each coefficient. The complete aerodynamic equations with definitions 
and explanations are provided in the appendix. 

1. Lift Coefficient is a function of: Mach No., altitude, angle-of- 

attack, leading-edge flap (LEF) deflection, trailing-edge flap (TEF) de- 
flection, horizontal tail deflection, speedbrake deflection, aileron 
deflection, pitch rate and angle-of-attack rate. 

2. Drag Coefficient is a function of: Mach No., angle-of-attack, 

LEF deflection, TEF deflection, horizontal tail deflection, aileron de- 
flection and speedbrake deflection. 

3. Pitching Moment Coefficient is a function of: same as lift co- 

efficient with the addition of rudder deflection. 

4. Yawing Moment Coefficient is a function of: Mach No., altitude, 

angle-of-attack, sideslip angle, LEF deflection, TEF deflection, differ- 
ential tail deflection, speedbrake deflection, rudder deflection, aileron 
deflection, roll rate and yaw rate. 

5. Rolling Moment is a function of: same as yawing moment with the 

addition of flaperon or differential TEF deflection. 

6. Side Force Coefficient is a function of: same as yawing moment. 
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C. AERODYNAMIC DATA 



Once the aerodynamic equations were obtained the next task was to 
obtain the value of each term in each equations for given flight condi- 
tions or, the aerodynamic data. This data, presented graphically [Ref. 4] 
was derived from wind tunnel testing but updated where possible by develop- 
mental flight test results. The data was given for low angle-of-attack 
and high angle-of-attack, considered to be forty degrees or higher. The 
distinction exists for the following reason: Above forty degrees angle- 

of-attack the leading-edge flaps are fixed to 34 degrees and the trail ing- 
edge flaps are undeflected. This is the configuration used in measuring 
the basic coefficients and no increments are added for leading or trailing- 

o 

edge flaps. Below 40 angle-of-attack the basic coefficients are measured 
at the zero flap deflections configuration and increments are added for 
leading and trailing edge flap deflections as necessary. 

Data was available for most of the flight envelope. In instances where 
no data was available, such as high angle-of-attack speedbrake data, the 
increments were set to zero. If for some increment data was not available 
throughout the desired ranges, judgment was made to determine the increment 
in one of three ways. 1) If the data reports noted that linear interpola- 
tion was possible, then the value was so obtained, 2) If it appeared that 
the increment was approaching to be zero, realistically it was made to go 
to zero or, 3) If no other indications existed, the increment was left con- 
stant throught the range. As an example, consider yawing moment increment 
due to speedbrake deflection. The data was presented for sideslip angles 
of positive two and ten degrees. The incremental changes were required 
over a sideslip angle range from negative twenty to positive twenty degrees. 
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The discrepancy was solved as follows. The increments were lineraly 
interpolated between zero and ten degrees and then held constant from 
ten to twenty degrees. The negative sideslip angle values were determined 
by using the negative of the positive sideslip angle values. These ad- 
justments to the actual aerodynamic data comprise a very small percentage 
of the data. They do not occur in any ciritical values of flight condi- 
tions and are determined realistically enough to have no adverse effect 
on the validity of the simulation. In contrast by covering the complete 
range of flight conditions, the aerodynamic buildup provides for a more 
realistic simulation. Once the aerodynamic data were obtained and evalu- 
ated, they had to be extracted from the graphical form to tabular form 
for computer entry. The values of flight conditions and surface deflec- 
tions for which data are tabulated is presented in the appendix. 

D. PROGRAMMING 

As originally envisioned, the aerodynamic buildup would be an integ- 
ral part of the main simulation CSMP program. This approach quickly ran 
into problems with handling functions of three and four variables, large 
quantities of numbers and sorting techniques. It was decided to program 
the aerodynamic buildup as a FORTRAN routine used as a subprogram in the 
flight simulation. The additional requirement of providing aircraft 
coefficients continuously throughout the flight envelope from aerodynamic 
data tabulated at specific intervals led to the use of a table look-up 
routine with interpolation functions, for intermediate flight conditions. 
The program procedure is exactly the same for each coefficient as 
follows : 
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1. A data file exists holding all the values of all the terms in 
the aerodynamic equation for given flight conditions. 

2. The program reads the data file and loads the data into program 
matrices. There then exists a separate data matrix for each term in the 
respective coefficients aerodynamic equation. 

3. These matrices are then printed out to display the data being 
used in the buildup. 

4. The program then makes calls to interpolation subroutines to de- 
termine the actual value of each term in the aerodyanmic equation for the 
existing flight conditions. 

5. The terms are then added appropriately to form the total coeffi- 
cient. 

An example of the program for lift coefficient is provided in the 
appendix. Each coefficient varies only in the terms of the aerodynamic 
equation. 

The aerodynamic buildup program is segmented into six major sections. 

Section One: Variable definition, explanations, declarations and 

program parameters. This section contains FORTRAN declaration and dimen- 
sion statements, program operation notes and control cards. The control 
cards provide the options for obtaining or deleting hardcopy output of 
the aerodynamic data, and the computed derivatives and coefficients. 
Additionally, the user can select to use test input flight conditions or 
inputs from another source, such as the main simulation program. 

Section Two: Aerodynamic data and constants. This section contains 

the read and write statements for each of the six coefficients to load 
the program data matrices and provide hardcopy output of the tabulated 
data , i f desi red. 
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Section Three: Test flight condition inputs. This section provides 

the operator the input of test flight conditions and control surface 
conditions. Standard day atmospheric tables incorporated in the program 
can also be selected if desired. This section can be totally deleted when 
the program is used as a subprogram to a larger simulation. 

Section Four: Aerodynamic buildup. This section contains for each 

coefficient the interpolation subroutine call statements to the data 
matrices to determine the value of the terms in the aerodynamic equa- 
tions. The total coefficients are actually determined in this section 
by summing the terms in the respective equations. 

Section Five: Output. This section contains the format statements 

for the formatted output of hardcopy data and results. 

Section Six: Subroutines. This section contains the interpolation 

subroutines used by the program. There are four subroutines, one each 
for functions of one, two, three or four variables. 
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V . SUMMARY 



There are no conclusions to draw for this report. A few comments 
can be made on the work that was done. The flight simulation project, 
initial project definition and planning was completed and the project 
is well underway. The initial concepts and requirements are still effec- 
tive, changed only for further clarification as work progresses. Ideas 
on programming are continually changing as problems are continually 
encountered and methods found to solve them. The final developed simu- 
lation will be completed in agreement with the ideas of this report. 

The aerodynamic buildup is complete. The program and data are on file 
at the Naval Postgraduate School, Monterey, CA. Point-of-contact is 
Dr. Marie D. Hewett, Department of Aeronautics (Code 67). The results 
of the aerodynamic buildup are verified as in agreement with tabulated 
and hand-calculated values. The programming though not extremely effi- 
cient, is simple, intelligible and flexible for use by various project 
members or even various projects. 
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APPENDIX A 



AERODYNAMIC EQUATIONS 
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REFERENCE T • E. FLAP DEFLECTION VALUES 
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APPENDIX C 



SAMPLE PROGRAM STATEMENTS 
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SECTICN U DEFINITIONS, EXPLANATIONS, DECLARATIONS AND 
PROGRAM PARAMETERS 
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PROGRAM CPERATI CK NGTES 
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(5) THIS AERODYNAMIC BUILDUP USES THE fGLLGHING SIGN 
CONVENTIONS FOR CONTROL SURFACE DEFLECTIONS 
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DIMEf'SIQI'/DECLARATION STATEMENT 
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CR1 (A, 18,11), CR2<A,2 ,18 , 11) , CR3 ) A, 2 , 18 , 11) , CR4t4,6,18), 
CRS(A»*)» CR6(A,2,11 ,18) , CR7 ( A ,3 , 11 , 1 fc ) , CR8<4,4), 

CR9 ( A , 5 ,1 8 ) , CR1 04 A, A ) , CR11(A,18), CR 12 i A ,4 , i 8) , CR13<4,4), 
CRIA(IE), CR1 5 (A, 2, 18), CR16<4,2), CR1 7 < A , A, 1 8 ) , 



CY1 (4 t J8* 1 1 1, CY2(4,2,18, 111 , CY3 i 4, 2 , 18 , ill » CY4<4, 6,181 
C Y5 ( 4 » 4 ) * CY6 ( 4 , 2 » 11 » 18) , CY7 i 4, 3 , 11 , 1 61 , CY8(4,4I, 
CYS(4,f,18lt CYl C( 4»4I * CYli(4,18l, CY12<4.18l t CY13(4,4) 
C Y 1 4 { 4 ,41 , CY 15( 4 , 41 
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SECTICN 2: AERODYNAMIC DATA AND CONSTANTS 
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MANfctVERING FLAP DEFLECTION < LtF ) - IVMF 
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LIFT COEFFICIENT - BASIC CONFIGURATION < CLBAS 



APPENDIX E 
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MAO NO. = 0.60 ALTO. = 0.0 ALFA = ALL 
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MACF NO. = 0.80 ALTO. = 20000.0 ALFA = ALL 
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